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Anisotropic D-dimensional polymerized phantom membranes are investigated within a nonper- 
turbative renormalization group (NPRG) framework. One focuses on the transition between a 
high-temperature, crumpled, phase and a low-temperature, tubular, phase where the membrane is 
flat along one direction and crumpled along the other ones. While the upper critical dimension 
_ -Due = 5/2 - is close to D = 2 the weak-coupling perturbative approach is qualitatively and 
quantitatively wrong. We show that our approach is free of the problems encountered within the 
perturbative framework and provides physically meaningful critical quantities. 
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Polymerized phantom membranes display a remark- 
able behavior as the temperature is varied [1, 2]. While 
crumpled at high temperature due to fluctuations, they 
exhibit a low-temperature flat phase with long-range ori- 
entational order [3-5] . The existence of an ordered phase 
in D = 2 despite Mcrmin- Wagner theorem originates in 
the existence of an anharmonic coupling between bending 
and elastic degrees of freedom or, equivalcntly, between 
out-of-plane capillary and in-plane phonon, modes [3]. 
This coupling induces a phonon-mcdiatcd long-range ef- 
fective interaction between the out-of-plane fluctuations 
that stabilizes a low-temperature flat phase for two - 
and even less - dimensional membranes [6-8]. Within a 
RG picture the coupling between the different modes is 
responsible for a stiffening of the bending rigidity con- 
stant at low-momentum - n(q) ~ q~v with r/ > 
called anomalous elasticity, that suppresses the destruc- 
tive fluctuations in the transverse direction to the mem- 
brane [3, 6, 8]. Note that, although strongly reduced, the 
out-of-plane fluctuations are not completely suppressed 
and the membrane still display ripples [9-11]. 

This phenomenon of long-range order induced by cou- 
pling between bending and elastic degrees of freedom has 
provided the theoretical grounds for the existence and 
stability of polymerized membranes and, incidentally, of 
the recently discovered graphene [12], the first example 
of truly two-dimensional membrane. It has also under- 
lined the importance of the nature of the membrane in- 
ternal structure regarding the kind of order displayed at 
long distance. As a remarkable illustration of this fact is, 
in marked contrast with the case of polymerized mem- 
branes, the absence of long-range orientational order for 
membranes deprived of fixed connectivity, and thus of 
elasticity modulus, whose prototypical example are liq- 
uid membranes. Past studies have considered a wide va- 
riety of alterations of the conventional crystalline order 
of polymerized membranes: presence of topological de- 
fects, like dislocations or disclinations, vacancies, local 
variations of connectivity, random impurities, . . . giving 



rise to novel and rich critical behaviors including new 
universality classes (see [1] and [13] for reviews). 

A fruitful modification of the translation- 
ally/rotationally invariant internal structure of the 
membrane is the inclusion of an in-plane anisotropy. In 
[14, 15] it has been shown that membranes with intrinsic 
one-directional anisotropy display, when the temperature 
is lowered, a tubular phase, fiat in one direction and 
crumpled in the other ones. This phase lies in a temper- 
ature range intermediate between a high-temperature 
regime, corresponding to a fully crumpled phase and a 
low-temperature regime, where the membrane adopts a 
flat conformation in all directions, see Fig.l. 
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FIG. 1: Flat, tubular and crumpled phases of anisotropic 
membranes as functions of the temperature. 

While the existence and stability of tubular phases 
have been numerically confirmed [16, 17], the crumpled- 
to-tubule transition remains numerically and experimen- 
tally unexplored. The main reason for this situation is 
the lack of accurate predictions for the critical quantities. 
This is a priori a rather surprising situation. Indeed, a 
striking consequence of the presence of anisotropy is that 
the upper critical dimension is lowered from D uc = 4, for 
isotropic membranes, to D uc = 5/2. From the proximity 
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of D uc with D = 2 one could expect the perturbative ap- 
proach performed in the parameter e = 5/2 — D to be in a 
good position to evaluate the critical quantities in D = 2. 
However, as shown explicitly by Radzihovsky and Toner, 
due to the fractional nature of the upper critical dimen- 
sion, the e-expansion is, on the contrary, " extremely unre- 
liable" and even " qualitatively wrong" [15]. This is clear 
in view of the exponent 77 describing the correlation be- 
tween the tangents of the membrane, which is found to 
be negative. A genuine negative value for this exponent 
would correspond to a downward rcnormalization - a de- 
creasing - of the bending rigidity constant, in contradic- 
tion with what is expected from physical grounds [15]. 
Moreover alternative perturbative methods used to inves- 
tigate the low-energy physics of membranes, e.g. the self 
consistent screening approximation (SCSA), that have 
successfully worked for isotropic membranes [18-20], ap- 
pear here hard to implement given the complexity of the 
field theoretical formulation of the model. 

In this article we investigate the crumpling-to-tubule 
phase transition of anisotropic membranes by means of 
a NPRG approach. This method [21] has been recently 
used to investigate the crumpling-to-flat phase transition 
and the flat phase of £>-dimensional isotropic membranes 
embedded in a d-dimensional space [22] (see also [23]). 
One of its important features is that it does not rely on 
the proximity of an upper or lower critical dimension. 
This is the reason why it has allowed to confidently in- 
vestigate the physics of membranes in the whole (D, d) 
plane [22]. Such a method is appropriate in the present 
context where the e-expansion displays a pathological be- 
havior. 

The NPRG approach relies on the use of a running 
effective action [21] (see [24-28] for reviews), Tfe[r], a 
functional of r = r(x) a d-dimensional external vector 
that describes the membrane in the embedding space, x 
being a set of internal D-dimcnsional coordinates label- 
ing a point of the membrane. The index k stands for a 
running scale that separates the high-momentum, with 
q > k, from the low-momentum ones, with q < k and 
Tk [r] represents a coarse grained free energy where only 
fluctuations with momenta q > k have been integrated 
out. The running of k towards the value k = thus 
corresponds to gradually integrate more and more low- 
momentum fluctuations. The fc-dependence, RG flow, of 
Tfc is provided by an exact evolution equation [21]: 

^ = X (r ™ + fl ' r '^r} (1) 

where t = lnfc/A, A being some microscopic, lattice, 
scale. The trace in (1) involves a Z?-dimensional momen- 
tum integral as well as a summation over vectorial in- 
dices. The function Rk(q) realizes the split between low- 
and high-momentum degrees of freedom. Several forms 
of Rk{q) will be considered in the following. Finally Fi ' 



is the second functional derivative of [r] with respect 
to r. 

The effective action r^fr] relevant to study polymer- 
ized membranes must be invariant under the group of 
translations, which implies that it only depends on pow- 
ers of the tangent vectors c^r. For anisotropic mem- 
branes rotational invariance within the membrane is ex- 
plicitly broken between one direction - named y - and 
the remaining D — 1 dimensions which are kept isotropic. 
The effective action thus writes: 

r*[r]= [d^xdy /^r) a +tx(#r) a 

1 1 (2) 

where Z y , u v , ij_ and Cy & r e running coupling constants: 
Z y is a bending rigidity, u y an elastic constant and t± a 
tension term in the Indirections (indexed by /i that runs 
from 1 to D — 1). This last term, that plays the role of a 
temperature in the ^-directions, is kept non- vanishing 
since these directions are not critical. Finally C, y , or 
rather the temperature along y, t y = u y C y , parametrizes 
the approach to criticality in the y-direction. Note that, 
up to change of coupling constants, this action is that 
used in [14, 15]. Let us consider the mean-field treat- 
ment of this model for u y > 0. When ( y = 0, the min- 
imum of Tfe is given by a configuration where d y r van- 
ishes which characterizes a crumpled, disordered, phase. 
When Q y > 0, this minimum is given by a configura- 
tion r(x) = ( y y e D where e D is a unit vector that spans 
the one-dimensional submanifold along which long-range 
order occurs. Thus, when the temperature t y is low- 
ered, the system undergoes a phase transition between 
a crumpled phase at high temperature, with Q y = 0, 
and a tubular phase at a low, intermediate, tempera- 
ture, with Q y 7^ 0. When the temperature is further 
lowered, one recovers an isotropic flat phase along all the 
£)- directions, see Fig.(l). A specificity of the model is 
that at the crumpling-to-tubule transition power count- 
ing leads to q± oc q y , an anisotropic scaling which char- 
acterizes a Lifshitz-type behavior [29]. According to this 
scaling behavior the most relevant terms are those enter- 
ing in Eq.(2). Note that since we treat nonperturbatively 
this action we are not supposed to base our computation 
on power-counting arguments. We do this nevertheless 
since i) the proximity of D = 2 with the upper criti- 
cal dimension D uc = 5/2 leads to guess that the terms 
present in Eq.(2) play the major role ii) we wish to com- 
pare our approach with the previous perturbative ones 
[14, 15]. The assumption i) must, of course, be checked, 
which can be partly done by evaluating the sensitivity 
of the results with respect to modification of the cut-off 
function Rk{q). 

Let us now define the critical quantities. First, due to 
the anisotropy between the y-dircction and the remain- 
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ing D — l _L-directions one has the scale transformations: 
X _]_ — k x i and y = k z y' that define the anisotropic scal- 
ing exponent z. Under a RG transformation the field r 
acquires an anomalous dimension 77 such that Z y ~ k~ r ' z , 
while t± ~ fc _??i so the field r generically scales as: 
r = ) { 2( D - 1 - 3z + , ] z ) r' . One deduces the following rela- 
tion between z, rj and r/±: z = (2 — f?±)/(4 — ?/). Finally 
one defines the exponents z/ y and from the correla- 
tion lengths near criticality: £ y oc t y y and £j_ oc t~ v± 
with w y = zi/j. [15]. The flow equations for the cou- 
plings Z y , u y , t± and C, y are obtained using their defini- 
tions in terms of functional derivatives of the effective 
action and applying the RG Eq.(l). One defines the 
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dimensionless quantities: ( y = ky Z y 2 t ± 2 and 

-D-5 1 — D 

Uy = ky D ~ 5 Z y 2 t ± 2 u y with ky = k z . In terms of these 
quantities the RG equations write, with t — laky/ A: 

d t U = 



d t C y = -(2D-3-r,(D-3)/2)C 



HUy 



, . -2,D -2,D 

(d — 1) I 3-D q — 3 /q 3-D 



{5 -2D -77(5- D)/2)u y + 



(3) 



(£>-3)^ [(d-l)lt^; , + 9^f-^ 
while the anomalous dimension function r\ writes: 
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{D-3)u 2 y Cy - 108 Zp' 5-0 - 12(rf- llfs'-D p 



2,-D 



-540(^-5)^4^^ 

+ (D-5)(L>- 7) (-288^ C^o!^ +9M|j;^i 



+(d - 1) Mi'- 1 !) - 36 u y d N 6 q \o-d 



(4) 



, . . 7a , -D —ft . D . ^r-^a , D . . 

In Fqs.(3) and (4) l a b , N ab and M ab are the de- 
dimensionalized "threshold functions" (see [25, 26]) I"'® , 
■TV"'/ 3 and M" b D that are given by : 
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where = (7r/2) 
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[P(ly)] a [P(qy) 



(5) 



r[(3-D)/2],P(g„) = Z v gJ 



Rk(q y ), my = 4it, / £ y ! and where d/dt only acts on 
In Eq.(5) the function F[q y ] is given by 1, dP/dq 2 and 
(dP/dq 2 ) 2 for Z, JV and M respectively. These functions 
control the relative weight of the different modes: the - 
single - phonon mode of mass m y and the d—l capil- 
lary modes of zero mass, along the RG flow. Note, that 



these functions only depend on q y . Indeed enters 
only quadratically in the action (2) so that one can ex- 
actly perform the integration on the D — l ^-degrees of 
freedom in the /3 functions. 

Let us now discuss our results. First we find no renor- 
malization of tj_, i.e. rj± — 0, for any value of D, in agree- 
ment with the all-order perturbative result [15]. Due to 

the one- loop structure of the RG equation (1) we recover 

2 

the one-loop /3 function for u y and t y = u y ( y found in [15] 
using our equations (3) expanded around D = 5/2 — e: 



-2Cv+o(d+ 2 )*i 
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d t u y = -2 e u y + (d + 8) K 5/2 u 2 y 

with, at this order, a vanishing anomalous dimension rj* . 
Thanks again to the one-loop structure of (1), one gets 
the large-d results from the set of equations (3)- (4), with- 
out almost any more computation. In this limit, assum- 

2 

ing that u y ~ 0(1/ d) and ( y ~ 0(d) one finds, in agree- 
ment with this assumption, and using, e.g., the cut-off 
function Rk(q y ) — Z y (k 4 — qy)9(k 2 — q 2 ), a stable fixed 
point with coordinates: 



Ad (3-D) 
T (2D - 3) 



(5 - 2D) 



(6) 



Ad (3- D)(5- D)K D 



This fixed point exists for all values of D between 5/2 
and the lower critical dimension D = 3/2 (see [15]). 
The corresponding critical exponents are given by: v = 
1/(2D - 3) + 0(1/ d) and r\ = 0(1/ d). 

Going to finite values of d one finds a stable fixed point 
for any value of d in contrast to the isotropic case where 
it occurs only for d greater that some critical value [5] 
d cr (D) determined in [22]. One now specializes to the 
d = 3 case. One finds, for any dimension D < 5/2, a 
nontrivial anomalous dimension rj* from Eq.(4) at the 
fixed point. This exponent is found positive in any di- 
mension between D = 5/2 and D = 3/2, in agreement 
with what is expected. The critical exponents rj* and z/j_ 
in D = 2 are displayed in Figs.(2)-(3). They are repre- 
sented as functions of a real number A that parametrizes 
three families of cut-off functions among those that we 
have used: R\(q) = XR l k (q), i = 1,2,3 with Rl(q y ) = 
Zyk^M-qi/^), R 2 k (q y ) = Z y /(exp(q y /k i ) - 1) and 
Rl(q y ) = Z y {k 4 - ql)8(k 2 - q 2 y ). Variations of A al- 
low to investigate the cut-off dependence of the critical 
quantities and to optimize each cut-off function inside its 
family [30], i.e. to (try to) find stationary values of the 
critical quantities. For each exponent one succeeds to 
find a - single - stationary point for closed values of A, 
see Figs.(2)-(3). For u± one can extract from these con- 
siderations the value u±_ = 1.213(8) in agreement with 
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that found in [15]: v± ~ 1.227. For the anomalous di- 
mension one finds ry* ~ 0.358(4) which differs consider- 
ably from the value found in [15]: r/ ~ —0.0015. Finally 
one deduces from these results and the scaling relations: 
z = 0.5490(6) and v y = 0.665(5) to be compared to 
v y ~ 0.614. Note that variations of the critical quan- 
tities appear to be very smooth with both variations i) 
of the parameter A inside a family of cut-off function and 
ii) of the family itself, see Figs.(2)-(3). This large insen- 
sitivity with respect to the cut-off functions constitutes a 
strong indication of trustability of our results, in agree- 
ment with the proximity of the upper critical dimension. 
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FIG. 2: The exponent v± as function of the parameter A. In 
solid line R\(q y ), in dashed line R\{q y ) and in dot-dashed line 
R%(q y ). The dot on each curve corresponds to a stationary 
value of the exponent. 
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FIG. 3: The exponent 77 as function of the parameter A. The 
conventions are the same as in Fig.(2). 

We have shown that the NPRG allows to overcome the 
main difficulty that plague the perturbative, weak cou- 
pling, expansion for anisotropic membranes. Our pre- 
dictions can be easily tested through numerical investi- 
gations. This would validate the NPRG approach as an 
efficient alternative to the perturbative one, in particular 
in the perspective of probing the properties of membrane- 
like systems with various orders and geometries. 
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